Large Seebeck effect in the electron-doped FeAs 2 driven by quasi one dimensional 

pudding-mold type band 
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We investigate the thermoelectric propeties of the electron-doped FeAs2 both experimentally and 
theoretically. Electrons are doped by partially substituting Se for As, which leads to a metallic 
behavior in the resistivity. A Seebeck coefficient of about —200 pV /K is reached at 300 K for 1% 
doping, and about —120 /xV/K even at 5% doping. The origin of this large Seebeck coefficient despite 
the metallic conductivity is analyzed from a band structure point of view. The first-principles band 
calculation reveals the presence of a pudding-mold type band just above the band gap, somewhat 
similar to NaaiCoCh, but with a quasi-one-dimensional nature. We calculate the Seebeck coefficient 
using a tightbinding model that correctly reproduces this band structure, and this gives results close 
to the experimental observations. The origin of this peculiar band shape is also discussed. 

PACS numbers: 72.15.Jf, 63.20.dk, 71.20.-b 
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I. INTRODUCTION 



Searching materials with good thermoelectric prop- 
erties has been an issue of great interest both from 
the viewpoint of fundamental physics and application 
purposes. The efficiency of a thermoelectric material 
is often evaluated by the dimensionless figure of merit 
ZT = S 2 T/(pn), where S is the Seebeck coefficient, p 
the electrical resistivity, k the thermal conductivity, and 
T being the temperature^. S 2 / p in the ZT formula is 
called the power factor, and in order to have this quan- 
tity large, both large Seebeck coefficient and small re- 
sistivity are required. However, it is usually difficult to 
have large Seebeck effect in metallic materials, so good 
thermoelectric materials are often found in narrow gap 
semiconductors. 

The discovery of the large Seebeck effect in Na^CoO^ 
has provided a new avenue for searching good thermo- 
electric properties in materials with metallic resistivity 
with significantly amount of carriers. In fact, such a 
coexistence of metallic resistivity and large Seebeck co- 
efficient have also been observed in e.g., CuRhO^ 1 ^, 
LiRh 2 OA LaRhOA SrTiO^, and quite recently, in 
PtSb^- Theoretical studies on the origin of the large 
Seebeck effect in these materials have also been exten- 
sively performed^ - — . In particular, Kuroki and Arita 
pointed out in Ref. HH that a peculiar band shape re- 
ferred to as the "pudding-mold" type, which consists of 
a dispersive portion and a flat portion, can in general be 
favorable for the coexistence of large Seebeck coefficient 
and small resistivity. The quasi-two-dimensional pudding 
mold band has been proposed as an origin of the large 
power factor observed in Na^CoC^ and other materials. 
On the other hand, Mori et al. have recently studied the 
origin of the large Seebeck effect in PtSb2, and found a 
three dimensional "corrugated" flat band that gives rise 
to many Fermi surface pockets through the entire Bril- 



louin zone when holes are doped.— This also gives rise 
to a coexisting large Seebeck coefficient and metallic re- 
sistivity in general. 

In the present study, we focus on FeAs2 in an attempt 
to look for yet another material that exhibits large See- 
beck effect despite the metallic resistivity. FeAs2 crystal- 
lizes in an orthorhombic Marcasite structure with Pnnm 
space group (No. 58). The structure is characterized by 
the edge-sharing FeAs6 octahedra, which are connected 
along the crystallographic c axis (z axis) to form linear 
chains. The non-doped semiconductor FeAs2 has been 
investigated previously in Ref. [l^ in comparison with 
FeSb2, which exhibits a colossal thermoelectric effect^. 
A theoretical analysis has shown that the Seebeck coef- 
ficient of the non-doped FeAs2 can be understood quan- 
titatively using the band structure obtained from the 
first-principles calculation 2 ^. In the present study, we 
investigate the thermoelectric properties of the electron 
doped FeAs2 both experimentally and theoretically. By 
doping electrons by partially substituting Sc for As, the 
resistivity exhibits a metallic behavior, but nonetheless 
we observe a fairly large Seebeck coefficient at 300 K of 
about -200^V/K at 1% doping and -120 pV /K even 
at 5% doping. We also analyze the origin of this ex- 
perimental observation theoretically. It is found that the 
quasi-one-dimensional (qlD) version of the pudding-mold 
type band is present right above the band gap, and this 
peculiar band shape gives rise to the coexistence of the 
metallic resistivity and the large Seebeck coefficient. We 
also discuss the origin of this peculiar band shape us- 
ing a simplified model. With the present study, we now 
have examples of materials that possess one- (FeAs2), 
two- (Na x Co02, etc.), or three-dimensional (PtSb^) flat 
bands that give rise to the large Seebeck coefficient with 
the metallic resistivity. This reinforces the general effi- 
ciency of peculiar band shapes with partially flat dispcr- 
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II. EXPERIMENT 



A. Synthesis 



(a) 10" 



Polycrystalline samples of Fe(Asi_ 2; Se a ;)2 with nomi- 
nal Sc contents of x — 0.0, 0.01, 0.025, and 0.05 were 
synthesized by a solid-state reaction in two steps. First, 
stoichiometric amount of starting materials Fe (99.9%), 
As (99.9999%), and Se (99.9%) were mixed and ground. 
They were heated in an evacuated quartz tube at 500°C 
for 40 h and then at 700°C for 40 h. Second, the prod- 
uct was powdered, pressed into pellets, and sintered at 
800°C for 12 h. The obtained samples were character- 
ized by powder X-ray diffraction (XRD) and confirmed 
to be a single phase of Fe(As]_ x Se:c)2- Electrical resis- 
tivity p and Seebeck coefficient S were measured using a 
Physical Property Measurement System (PPMS, Quan- 
tum Design) in the temperature range from 2 to 300 K. 



B. Thermoelectric properties 

As shown in Fig. HJa), non-doped FcAs 2 exhibits a 
semiconducting behavior, consistent with the previous 
reports^ The Seebeck coefficient S of the non-doped 
FeAs2 exhibits a large maximum value of —0.65 mV/K 
at approximately 60 K, as shown in the inset of Fig.QJb), 
which is similar to the literature dataJ^ Then, both the 
magnitude and temperature dependence of resistivity p 
changed abruptly by substituting Se for As from semi- 
conducting to metallic. The p values on the order of 1 
mSlcm and the positive temperature slope of p suggest 
that a metallic sate is realized for Fc(Asi_ 2; Se a ;)2 (x = 
0.01, 0.025, and 0.05). The temperature dependence of 
S is also changed abruptly by Se doping. But even in 
the metallic state, the Seebeck coefficient reaches a fairly 
large value of about -200 pV/K at 300 K for the x = 0.01 
sample, as shown in Fig. [IJb). In this way, metallic p 
is compatible with a large S for Fe(Asi_a;Se x )2. This 
combination of the resistivity and the Seebeck coefficient 
gives a maximum power factor of 14 pW /cm. K 2 at ap- 
proximately 200 K for the x = 0.01 sample. 



III. THEORETICAL ANALYSIS 

A. Band structure 

We now turn to the theoretical analysis. We per- 
form first-principles band structure calculation using the 
Wien2k package 2 ^. The structural parameters adopted 
are given in Ref. in which the lattice constants are 
a = 5.2684,6 = 5.9631, and c = 2.9007 A, space group 
is Pnnm, and the atomic positions are (0,0,0) for Fe 
and (0.1760,0.3625,0) for As. In the present study, we 
construct tightbinding models in order to calculate the 
Seebeck coefficient. To obtain a model that correctly 
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FIG. 1: Temperature dependence of (a) electrical resis- 
tivity p(T) and Seebeck coefficient S(T) of polycrystalline 
Fe(Asi_ a; Se ;c )2 with x = 0.00, 0.01, 0.025, and 0.05. The 
inset shows Seebeck coefficient of non-doped FeAs2. 



reproduces the first-principles band structure, we con- 
struct maximally localized wannier functions (MLWFs) 
from the first-principles calculation result^ We focus on 
the electron doped case up to about the room tempera- 
ture, so that we only have to consider the band structure 
above the band gap. The tight-binding Hamiltonian is 
written as 



H 



ija \iv 



(1) 



where i,j denote the sites, and p, v the orbitals, and 
is the hopping parameter obtained from MLWFs , 
and c 1 , c are the creation-annihilation operators. Here we 
first consider six orbitals (three orbitals per Fe) to repro- 
duce the bands up to about 5eV above the Fermi energy. 
These orbitals have mainly 3d character, although they 
are hybridized with As Ap. 

The first-principles band structure is shown in FigUJa) 
with the tight binding bands superposed. The band 
structure just above the band gap has a nearly flat 
dispersion along U-R-T-Z (see (b) for the symmetrical 
points in the Brillouin zone), namely, within the k z = tt 
plane. In fact, the Wannier orbitals that are the origin 
of these bands are found to be elongated in the z di- 
rection (see FigfS]). We will call these orbitals the ll d z 2 
orbitals" for simplicity (the maximally localized Wannier 
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FIG. 2: (a) The band structure and the density of states of 
FeAs2- The solid black lines are the original first-principles 
band structure, while the red dashed lines are the bands of 
the tight-binding model constructed from MLWFs. (b) The 
Brillouin zone is shown, (c) A schematic image of the qlD 
pudding-mold type band of FeAs2. 



orbitals are obtained by projecting onto the d z 2 orbital), 
although they are actually more complicated due to the 
hybridization with the As 4p orbitals. The quasi-one- 
dimcnsionality of the bands within the k z = it plane is 
partially because of the anisotropy of these d z 2 orbitals, 
but to be more precise, there is another factor that re- 
stricts this quasi-one-dimensionality to near the k z = tt 
plane, which will be discussed later. Furthermore, these 
bands are somewhat flat also along Z-T (i.e., the k z direc- 
tion) up to some intermediate point, then bend sharply 
into dispersive portions. In Figl^c), we show a schematic 
figure of this band structure. It has a flat portion at 
the bottom with a qlD nature, so it is a qlD version of 
the pudding-mold type band introduced in reffLsl as the 
band structure of Na x Co02- This type of band structure 
is favorable for producing good thermoelectric properties 
when carriers are doped. The reason is because (i) even 
when a large amount of carriers are doped, the Fermi level 
sticks close to the band bottom due to the large density of 
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FIG. 3: The calculated Seebeck coefficient for x = 0.025 and 
x — 0.05 (the solid lines) compared with the experimental 
result (green squares and red circles). A calculation result 
with the band width reduced by 30% (green dotted line and 
red dash-dotted line). 



states, (ii) for such a position of the Fermi level, there is 
a large difference between electron and hole group veloc- 
ity, resulting in a large Seebeck coefficient, and (iii) small 
resistivity can be achieved owing to the large amount of 
doped carriers and the dispersive portion of the band 
near the Fermi level. The origin of such a peculiar band 
structure in FeAs 2 will be discussed later. 



B. Seebeck coefficient 

Let us now move on to the Seebeck coefficient calcula- 
tion. We first briefly introduce the method to calculate 
the Seebeck coefficient and the power factorji^ The See- 
beck coefficient S and the resistivity p are calculated with 
Boltzmann's equation using the tight-binding Hamilto- 
nian. They are given as, 



S 

P = 



J_Ki 

eT Kq ' 



(2) 
(3) 



where e is the elementary electrical charge(e < 0), T is 
the temperature. Kq and Ki are given as 



K„ = £V(kKv£ 



dm 

de 



(k) (4- M ) n , (4) 



where k is wave vector, r is quasi-particle lifetime, e£ 
is the v-th. energy eigenvalue at k, v£ = (l/h)(ds^/dk) 
is the group velocity tensor, / is the Fermi distribution 
function, and p is the chemical potential at T. We sim- 
ply denote {K n ) zz =_K n , S zz = S = (l/eT)(X /^i), 
and p zz — p = e~ 2 K . We will approximate the quasi- 
particle life time as an (undetermined) constant in the 
present study, so that it cancels out in the Seebeck coef- 
ficient calculation, while not in the resistivity. Thus the 
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FIG. 4: (a)The calculated Seebeck coefficient and (b) the 
normalized power factor at 300K against the doping level x. 



absolute value of the power factor S 2 /pis not determined 
and will be normalized by a certain value. 

The calculated Seebeck coefficient as a function of the 
temperature is shown in Figl3] The results at T — 300K 
give -100^V/K at x = 0.025, and -77/iV/K at x = 0.05. 
The theoretical result is about 20% smaller than the ex- 
perimental values at x = 0.025 and T = 300K. This 
discrepancy from the experiment may be due to the con- 
stant lifetime approximation, or the band width renor- 
malization due to the electron correlation effects that are 
not considered in the first-principles band calculation. 
As for the latter possibility, it has been known from the 
comparison between the band calculation and the angle- 
resolved photoemission studies that the bandwidth of the 
3d electron materials can be reduced by as much as 60% 
, and taking this effect into account in the Seebeck co- 
efficient calculation reproduces quantitatively the exper- 
imental results of Na^CoO^. The Seebeck coefficient 
is underestimated also in other cases of 3d electron sys- 
tems such as SrTiO;?i£, due to the band width renor- 
malization effect. If we assume that the band width is 
reduced by 30% from the first-principles result, the cal- 
culated Seebeck coefficient is in good agreement with the 
experimental result as shown in FigJIJJ dotted line). The 
band width renormalization effect is not so strong as in 
Na^CoC^, indicating that the electron correlation effects 
are not so large in this material. 

In FigHJ we show the doping level dependence of the 
Seebeck coefficient and the normalized power factor (nor- 
malized by the maximum value) at T = 300K. This re- 
sult shows that the Seebeck coefficient keeps relatively 
large values of about — 50/xV/K even at a (hypothcti- 




FIG. 5: (a) A simplified tightbinding model of FeAs2 and (b) 
crystal structure of FeAs2- The band structure in the unfolded 
Brillouin zone is shown in the right panel with (c)ti z = — 1, 
(d)ti z = —1 and t2z = —0.4, and (e)ti z = —1, ti z = —0.4 and 
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-0.1. 



cally) large doping rate x ~ 0.5, so that (assuming a con- 
stant quasiparticle life time) the power factor is roughly 
constant. This peculiar doping level dependence of the 
power factor is a consequence of the pudding-mold type 
band. 



C. Origin of the pudding- mold type band 

As seen from our theoretical analysis, the large Seebeck 
coefficient in the electron-doped FeAs2 is a consequence 
of the qlD pudding-mold type band. In this section, we 
focus on the origin of this peculiar band structure. For 
simplicity, we first consider a unit cell that contains only 
one Fe site per unit cell. This corresponds to expanding 
(unfolding) the Brillouin zone. By taking this unit cell, 
we consider a simplified tight-binding model, where only 
the hoppings t z \ and t Z 2 along z direction, and t xyz to 
the (1/2,1/2,1/2) position are considered (FigJSJa), see 
FigGJb) for the correspondence with the actual lattice 
structure). The hoppings in the x and y directions can 
be neglected for the band originating from the d z i orbital. 
In this model, the energy is given as 

E(k)= - 2t lz cos(k z ) - 2t 22 cos(2fc z ) 

— 8t xyz cos(k x /2)cos(k y /2)cos(k z /2). (5) 
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FIG. 6: (a) The two orbital tight binding model constructed 
from MLWFs and (b) the band structure of Fig[5jc) shown in 
the folded Brillouin zone. 



and the distant hoppings are also contained. The band 
structure of the two orbital model is shown in Fig(5la) 
compared with the band structure of Fig[5jc) refolded 
into the original Brillouin zone (Figj6jb)). The band 
structure of the simplified model is roughly similar to 
that of the two orbital model especially near the band 
bottom, so that the essence of the origin of this pecu- 
liar band structure can be understood by the simplified 
model. To be precise, the magnitude of t 2z and t xyz is 
different between the two models because the simplified 
model considers only three hoppings, while the more re- 
alistic two orbital model takes into account the distant 
hoppings as well. The bottom line here is that the origin 
of the qlD pudding-mold type band is the overlapping 
feature of the d z i orbitals, where a relatively large sec- 
ond neighbor hopping is present. 

IV. CONCLUSION 



When k z ~ n, the third term in the right hand side van- 
ishes due to cos(7r/2) = 0, and the band dispersion loses 
k x and k y dependence. This is the explanation of the qlD 
nature of the band structure in the k z = ir plane. We now 
focus on the band structure along the k z direction. When 
we consider only t\ z , the band shape is free-electron-likc 
at the bottom of the band(Figj5jc)). Adding the distant 
hopping t 2z = —0.4, the shape of the band becomes flat 
near the band edge around k z — it (corresponding to 
the Z point), so that the qlD pudding- mold- type band 
is formed (FigJ5Jd)). Further considering t xyz = —0.1, 
the band becomes asymmetric with respect to k z = tt 

(FigEte)). 

To compare this simple band structure to a more re- 
alistic one, we construct, again exploiting the MLWF's, 
a two orbital model that considers only the d z 2 orbital 
in the original unit cell that contains two Fe. Although 
this model docs not perfectly reproduce the original band 
structure as in the six orbital model presented in the 
above, it still captures the essential features. In this 



In the present work, we have studied the thermoelec- 
tric properties of the electron-doped FcAs2 . A large See- 
beck coefficient is observed despite the metallic behavior 
of the resistivity. First-principles band calculation re- 
veals the presence of a qlD pudding-mold type band just 
above the band gap, which mainly originates from the Fe 
3dz2 orbital, orbital. Using a tightbinding model that 
correctly reproduces this band structure, we have calcu- 
lated the Seebeck coefficient. The results are close to 
the experimental observation, and the electron correla- 
tion effects are found to be small. The qlD nature of 
the pudding-mold type band comes from the combina- 
tion of the d z 2 orbital and the lattice structure. Given 
the present study, we now have examples of materials 
that possess one-, two-, or three-dimensional flat bands 
that give rise to the large Seebeck coefficient in a metallic 
system. Therefore, the present findings in the electron- 
doped FeAs2 reinforces the general efficiency of peculiar 
band shapes with partially flat dispersion. 
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